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Abstract 


In this paper we study a real scalar field as a possible candidate to explain the dark matter in 
the universe. In the context of a free scalar field with quadratic potential, we have used Union 2.1 
SN la observational data jointly with a Planck prior over the dark matter density parameter to set 
a lower limit on the dark matter mass as m > eV {c = h = 1). For the recent value of the 

Hubble constant indicated by the Hubble Space Telescope, namely Hq = 73 ± 1.8 km s“^Mpc“^, 
this leads to m > 1.56 x 10“^^ eV at 99.7% c.l. Such value is much smaller than m rsj 10-22 eV 
previously estimated for some models. Nevertheless, it is still in agreement with them once we have 
not found evidences for a upper limit on the scalar field dark matter mass from SN la analysis. In 
practice, it confirms free real scalar field as a viable candidate for dark matter in agreement with 
previous studies in the context of density perturbations, which include scalar field self interaction. 
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I. INTRODUCTION 


Nowadays the most accepted model in cosmology is known as ACDM model, where CDM 
stands for Cold Dark Matter and A is the cosmological constant term, the latter being the 
main candidate to explain the current phase of acceleration of the universe and the former 
having a central role for the structure formation in the standard cosmology. The pressureless 
matter represents about 30% of the total material content of the universe, where about 5% 
is baryonic matter and 25% is nonbaryonic dark matter (DM). The cosmological constant 
stands for the 70% remaining part, sometimes also associated to a dark energy (DE) exotic 
fluid (see [l| for a review). 

The idea of an accelerating universe is indicated by type la Supernovae observations [2|, 
in agreement with ACDM, but the nature and origin of the dark matter is still a mystery 
(see^, 1^ for a review and for a ten-point test that a new particle has to pass in order to 
be considered a viable DM candidate). 

The hrst candidate to DM as a scalar held is the axion, one of the solutions for the 


Charge-Parity problem in QCD [71,18|. The axion is essentially a scalar held with mass of 
about 10“® eV, which has its origin at 10“^° seconds after the big bang. This candidate is 
until now one of the most accepted candidates to DM particles. Among others candidates 
for DM particles in the universe, the so called scalar held dark matter (SFDM) model 0-Q 
is one of the most studied models in quantum held theory, and its applicability in cosmology 
to explain the diherent processes in the evolution of the universe has been investigated in the 
last decades (see 15| for a review and references therein). In 1^ it is proposed a cosmological 


scalar held harmonic oscillator model, in agreement to ACDM model. Scalar helds are also 
used to drive the inhationary phase of the universe [l3, [l^. One of the hrst applications 
of a complex scalar held for structure formation of the universe was given by Press and 


Madsen 19|, l20|]. Most recently, in 2l|] it was studied structure formation by assuming that 


dark matter can be described by a real scalar held. It was also investigated the symmetry 
breaking and possibly a phase transition of this scalar held in the early Universe. At low 
temperatures, scalar perturbations of SFDM leads to formation of gravitational structures. 


In 22|] it has been shown that such model agree with rotation curves of dwarf galaxies 
and small and large low surface brightness (LSB) galaxies. The model was also extended 
to include the dark matter temperature and hnite temperature corrections up to one-loop 
in the perturbative regime. Gravitational constraints imposed to dark matter halos in the 
context of hnite temperature SFDM has been investigated in 23|]. In 2^ it was studied the 
Sextans dwarf spheroidal galaxy, embedded into a scalar held dark matter halo. In 1^ the 
quantum creation of scalar particles was studied in the context of an alternative accelerating 
model of the universe. 


It is important to notice that some models 13|, llJ, 12^ points to a mass of the scalar held 


of about 10 eV, while a recent work points to 10 eV. As we shall see, our model 
is in agreement to the last one. 
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In this work we study a real scalar field evolving in a Friedmann-Robertson-Walker back¬ 
ground as a possible candidate to explain the dark matter in the universe. In section II, 
we develop the dynamics of the model and analyse the stability of the resulting system of 
equations. In Section III, we make an alternative derivation based on a change of variables, 
more suitable for numerical integration. In section IV, we constrain the free parameters of 
the model by using observational SN la data. We conclude in Section V and deduce the 
scalar held motion equation in the Appendix. 


II. THE DYNAMICS OF THE MODEL 


We are interested in the study of the action of the form 




^grav _|_ ^ 


mat 


( 1 ) 


where 

I + 2A) (2) 

is the standard Einstein-Hilbert action for the gravitational held, R stands for the Ricci 
curvature scalar, A is the cosmological constant parameter and G is Newton’s gravitational 
constant. 

The second term in ([T]) stands for the matter content of the universe, namely, 

Smat ^ j + j , ( 3 ) 

which is composed of the Lagrangian density of the standard matter contributions (radiation, 
baryons, neutrinos etc.), in addition to the Lagrangian density of a real massive scalar 
held 0 minimally coupled to gravity, given by 


c* = - V(4 >), (4) 

where V (0) is the potential. 

Variation of the action S with respect to 0 leads to the equations of motion for the real 
scalar held. The equations of motion for the gravitational held are obtained by varying the 
action with respect to the metric g^'^, which leads to 


Rfiu 2^9^i/R ^'xGT^y , 

where the energy-momentum tensor of the matter helds, is dehned as 


T 


2 55 ™“* 


rjism 


-f T'^ 


(5) 

( 6 ) 
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Here we are assuming that the standard matter term in the above equation can be ap¬ 
proximated by a perfect fluid characterized by its energy density p*™', pressure and 
four-velocity Therefore, we have 

T;™ = {r + . (7) 


On the other hand, the scalar field part of the energy-momentum tensor is 

7^ = - V (0)] , (8) 

which can be rewritten in the form of a perfect fluid ([7]) using the definitions 

p^=^d°‘(/)da(/> + V((/>), = ^d°‘(/)da(/> - V((/)), (9) 

In particular, assuming that the field is homogeneous (9,0 = 0), we have 


= + v{4>), = ^ 0 ^ -'' 7 ( 0 ) • ( 10 ) 

We restrict our treatment to a universe composed of a single baryonic component with 
energy density pb, a cosmological constant term with energy density pA, and the scalar field 
dark matter component p^j,. The equations of state of those components are respectively 
Pb = 0, Pa = -PA, and p^ = oj^p^, where 


n _ W ~ ^ 

P0 “ 102 + v{ct)) ■ 


The Friedmann equations obtained from this approach are (in units c = h = 1). 


with the Friedmann constraint 



(12) 

■ dV 

0 + 3770+ — = 0, 
d<p 

(13) 

Pb + 3 77pfe = 0, 

(14) 

Pa = 0, 

(15) 

^2 

77^ = — (p0 + Pb + Pa) , 

(16) 


where = SvrG and H = a/a is the Hubble parameter. 

In order to study the system of equations ffT2l) - ffT5ll , we define the following dimensionless 
variables (accordingly to the variable transformations adopted in to a similar set of 
equations): 


K 0 


K 




a/6 77’ 


u = 


W 


\/3 H 


h = ^ 

■ a/3 77 ’ 


_ K a/^ 

a/3 77 • 


(17) 
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Consequently, the Friedmann constraint ffT 6 |l reduces to 

2 I 2 I 7 2 I 72 -1 

X + U +0 +/ =1. 


( 18 ) 


Specializing to a quadratic scalar potential C(0) = where m is the mass of the 

scalar held, the background evolution of the universe is given by the equations as 


x' = —3x — su + -Ilx, 

3 

u' = sx + 

h' = |(n-i)(>, 
r = |ni. 
s' = 5ns. 

where a prime denotes derivatives with respect to the e-folding number = In a. 


(19) 

( 20 ) 
( 21 ) 
( 22 ) 
(23) 


s = m/H , 


(24) 


and 


3 H 3 2 72\ 

2" " “in = 2^“^ + * 5 


(25) 


In order to study the stability of the system ffT^ - fl23l) from a dynamical system approach 


26l |. we dehne the vector v = (x, u, b, I, s) and consider a linear perturbation of the form 


V ^ V + Sv. Thus, the linearized system reduces to Sv' = MSv, where M is given by 


-3 9x2 1^2 

—s 3bx 

6xu + s 

|62 + 3a;2 3bu 

6 x 6 

0 |62 -I- 3x‘ 

6x1 

0 

3bl 

6xs 

0 

3bs 


0 \ 
0 X 

0 0 

|62 + 3x2 Q 

0 §62 + 3x2 y 


(26) 


and represents the Jacobian of v'. 

The equilibrium points or hxed points {xc, Uc, be, h, Sc} of the phase space are: 

(I) {±1, 0, 0, 0, 0}; 

(II) {0, 0, ±1, 0, 0}; 

(III) {0, ±VT^, 0, /, 0}; 

(IV) { 0 , 0 , 0 , 1 , 4 . 

According to the Friedmann constraint flT 8 |l . the case (I) represents a totally kinetic energy 
scalar dominated universe, the case (II) characterizes a totally baryonic dominated universe, 
the case (III) represents a mixture of scalar potential energy contribution and cosmological 
constant and the case (IV) a totally cosmological constant dominated universe. The case 
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(IV) is the only one that admits s > 0, leading to a mass m > 0 for the dark matter particles 
from flM)) . The eigenvalues of the matrix M at the hxed points above are: 

(I) { 6 , |, 3, 3, 3}; 

(JJ) 3 3 S'). 

(III) {0, 0, 0, -3, -§}; 

(IV) {0, 0, -| + -| - -§}. 

As it is well known from the dynamical system approach [26|, when the real parts of the 
eigenvalues are positive the hxed point is unstable, and when the real parts are negative the 
hxed point is stable. When positive and negative real terms are present, the hxed point is a 
saddle point. Finally, when some of the eigenvalues are null, nothing can be said about the 
stability, and more accurate methods should be applied, as the center manifold theory. 

For the above set of eigenvalues, the most interesting one for us is the case (IV), where 
the stability is not completely discarded due to the null values of the two hrst eigenvalues. 
Beside that such case admits s > 0, leading to a non-null mass to the dark matter particle. 
Although such case corresponds to a completely cosmological constant dominated universe 
at the hxed point, it is possible to estimate the mass as m Ho 10 eV if we assume 
s ~ 1. The case (III) is also interesting, since it also has some of the eigenvalues negative 
and others null, but it corresponds to s = 0, leading to m = 0. This case is also interesting 
since it admits the possibility of coexistence of cosmological constant and a potential part 
of <p. 

The limit of large mass can also be studied. From flTT]) . we have —1 when m —)• oo 

or equivalently V(0) >> 0^/2, which corresponds to x « u. Taking the limit a: = 0, the 
hxed point of interest of the dynamical system is {0, 0, 0, 1, s}, indicating that the large 
mass limit corresponds to the ACDM model, with I = 1 and s —)■ oo. 

In order to study the above system quantitatively we should try to integrate the dynamical 
system (HlD-dSSl). Nonetheless, we will do an alternative treatment in the next section, aiming 
the numerical integration. 


III. ALTERNATIVE DERIVATION - CHANGE OF VARIABLES 


The system of diherential equations derived in last section is useful for dynamical system 
analysis, as we had discussed in the last section. In order to constrain the SFDM mass 
we should solve these equations numerically. However, we must stress that this choice of 
variables may generate high frequency solutions in the variables x and m, as already presented 
in Ref. 13|] and shown in Fig. [H For huge values of the mass m, these high frequency 
solutions require high computational efforts (as an illustration, one single evaluation for 
m ~ takes more than 24 hours with an Intel i5 processor), being desirable to avoid 

such parametrization. 

Given that the behaviour of and Vtp = x^ + are smooth even for large values of 
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FIG. 1: Panel a) Evolution of variables {x,u,b,l) in function of redshift for /r = 10^ and uq = 0. 
Panel b). Evolution of variable s in function of redshift for /r = 10^ and uq = 0. 

the mass m, a reasonable choice of variables is r = such that and the 

“angle” 6 = tan“^ Thus, the new variables are: 


X = r cos 9 
y = r sin 6 . 


(27) 

(28) 


In some sense we are changing from cartesian to polar coordinates. The resulting system of 
equations is now: 


3 

—3r cos^( 6 ') + 2 ^^ 

(29) 

s + ^sin 2 ( 0 ). 

(30) 

1 — 1 

1 

CO 1 CM 

(31) 

|n;. 

(32) 


(33) 


where 11 = 2r^cos^(0) + 6 ^. As can be seen, the equations are now even simpler on r and 
6, and the system is invariant under the transformation 6* —)■ 6*i + tt, in such a way that it 
is periodic over 6 with a period tt. Thus, we seek to solve this system of equations over the 
vector of variables r = (r, 6, b, I, s) with the corresponding initial conditions at = 0 : 


fh — 9o, \/f2fto, y/f — hj 


(34) 


where Oq becomes a free parameter, Oq e |], we have used the normalization condition 
Qrn + = 1 and we have dehned the parametrization of the mass: 

m 


y = 


Hn 


(35) 
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FIG. 2; Panel a) Evolution of variables (r, 6, 1) in function of redshift for /r = 10^ and 9q = 0. Panel 
b). Evolution of variables (s,0) in function of redshift for fj, = lO”^ and 9o = 0. We have plotted 
9 + ^ here in order for both variables stand on same axis scale. 

It is also interesting to notice that in these new variables, the oo^ EOS assnmes a nicely 
simple expression: 

= cos(26*) (36) 

In Figure El we may see that none of the “polar” variables present oscillations for higher 
mass values (/i = 10^ on Fig. E])- 


IV. CONSTRAINING DARK MATTER MASS 


In order to obtain limits to the dark matter mass, we have to constrain the set of free 
parameters of the SFDM model, namely, s = {Q^o,9o,flbo, Hq). Hence, even imposing 
spatial flatness, as predicted by inflation and indicated by CMB observations, we end up 
with 5 free parameters. Among these parameters, Ofeo is the most well constrained to date, 
with compatible limits from BBN and CMB. We therefore choose to £x it at the value given 
by Planck and WMAP: Om = 0.049 27|, reducing the number of free parameters to 4. 

As a hrst attempt to constrain these free parameters we have considered the measurement 
of the Hubble parameter H{z) in different redshifts. These kind of observational data are 
quite reliable because in general such observational data are independent of the background 
cosmological model, just relying on astrophysical assumptions. However, even using the 
current most complete compilation of H{z) data, with 34 measurements 28|], it was not 
enough to constrain the 4 free parameters, mainly because of the oscillatory behaviour of 
the model along with the relatively small number of degrees of freedom [u = n — p = 30). 

Next, we considered a SNe la data sample, which, although being more dependent on 
the hducial cosmological model, it consists of a large data sample and has passed recently 
through more rehned and model independent methods of light curve htting 29|. 












The parameters dependent distance modnlns for a snpernova at a redshift 2 ; can be com- 
pnted throngh the expression 


I^sn(~-\p) = msN - Msk = 5 logdi + 25, 


(37) 


where msN and Msn are respectively the apparent and absolnte snpernova magnitndes, 
p = fi, Hq) is the set of free parameters of the model and di is the Inminosity 

distance in nnits of Megaparsecs. 

Since we have no analytic expression for H{z), it is necessary to define di throngh a 
differential eqnation. The Inminosity distance di can be written in terms of a dimensionless 
comoving distance D by: 

d^ = (l + z)^D. (38) 

c 

On its tnrn, the como ving distance can be related to H{z), for a spatially flat Universe, 
by the following relation 30|: 

dB 1 

(39) 


dP 

dz 


E{zy 


where E{z) = ^ 7 ^- Changing to the independent variable N and given that E = we 
arrive to: 

p-N 

D'{N) = (40) 

/i 

which can be seen as a sixth eqnation to be solved simnltaneously with the system ([29])-([33]) 
with the initial condition for the comoving distance, T*(0) = 0. In order to constrain the free 
parameters of the model we considered the Union 2.1 SN la dataset from Snznki et ah 29]. 
The best-£t set of parameters p was estimated from a statistics with 

[fJ-SNiZilp) — PsN,o,if 


XSN — 


2=1 


a? 


(41) 


where pswi^ilp) is given by fl37|l . psN,o,i is the corrected distance modnlns for a given SNe la 
at Zi being Uj its corresponding individnal nncertainty and n = 580 for the Union 2.1 data 
compilation. 

As nsnal on this analysis, we marginalize over the EIq dependence by rewriting the distance 
modulus: 

^-SNiz) = blogDiiz) + M* (42) 

where Dl = {1 + z)D is dimensionless luminosity distance and M* = 25-1-5 log comprises 
all the dependence over Hq. Then, we marginalize the likelihood over M*: 


^+00 


C{Vt^o,do,pi) = / AAexp 


2^ (Af*, 12,^0) ^0) h") 


dM^ 


(43) 


where AA is a normalization constant. The corresponding = —2 In is given by: 


( 44 ) 
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FIG. 3; Confidence contours of SFDM from Union 2.1 SN data. The regions correspond to 68.3%, 
95.4% and 99.7% c.l. (green, red and blue regions, respectively). 


^\og[DL(zi)\-^lo,i 




b\og[D L{zi)]- 


Where 71 = = Er=i 

Since we are mainly interested on the constraints over the dark matter mass, we have 
marginalized C over the parameter 9q. The result of this analysis can be seen in Fig. |3l 
where we have plotted the statistical conhdence contours on the plane logj^o h ‘ 

As one may notice, there is a strong degeneracy for in the range of low DM mass 
(10“^ ^ /i ^ 10). For higher masses, fi ^ 10^, the degeneracy changes to prevent a mass 
determination. 

In this case, restricting to the quite wide range 0.04 < < 0.6 (flat prior on this 

interval), we have found a xLin = 560.854, which corresponds to xl = 0.974. We have found 
the best £t = 0.22^Q;^5^Q;^7^Q;^g at = 2.30, 6.17 and 11.83, respectively. The best 


--0.15-0.17-0.18 


£t for n was 19.95 and we may infer lower limits to /i as 0.94, 0.51 and 0.13 at the same 
levels. 

In order to alleviate the degeneracy over we use a prior over based on Planck 
results in order to constrain the mass m, or, at least, give it an inferior limit. As Planck 
(2013) [27| yields the limit fldmo = 0.265 ± 0.011, at la, we use a Gaussian prior of = 
0.265 ± 0.022, which is approximately 2a of the Planck analysis. The result is shown on 
Figure HI 

As one may see, is now better constrained in the interval 0.19 < < 0.34, at 

3a. We have found, from this analysis, Xmin = 561.897, xl = 0.976. The best £t is 
= 0-252+Q;Q27lao45^o;o6i- The degeneracy remains, however, over the dark matter mass, 
whose likelihood now presents two maxima, one for low mass and another for high mass. 
The best £t is /i = 19.68, with inferior limits given by 1.12, 0.64 and 0.12 at = 2.30, 
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FIG. 4: Confidence contours of SFDM from Union 2.1 SN data + Planck prior over U^o- 
regions correspond to 68.3%, 95.4% and 99.7% c.l. (green, red and blue regions, respectively). 


The 


6.17 and 11.83, respectively. 

The lowest limit on this context, /r = 0.12, corresponds to a dark matter mass m = 
0.12i7(C^, which corresponds to m = (1.557 ± 0.038) x 10“^^ eV for a Hubble constant of 
Hq = 73 ±1.8 km s“^Mpc“^, as indicated by the most recent analysis from the Hubble Space 
Telescope [31 1. 


V. CONCLUDING REMARKS 


In this work, we investigated the hypothesis of the dark matter as a scalar field, where 
we have analysed the dynamics of a quadratic potential for a single parameter free scalar 
held. Considering the SN la observational data, we have found a very low inferior limit to 
the dark matter mass, namely, m > 1.56 x 10“^^ eV. Although such lower limit is much 
smaller than the constraints obtained in 00,0 (of about 10 eV), we have shown 
that any dark matter mass greater than this one is also compatible with SN la observations. 
Our constraints are also much similar to the value recently obtained in 0 (m ~ 10 eV), 
which opens the question if the mass of the scalar held dark matter could be so small. 

A very small mass may be in conhict with the structure formation of the universe, since it 
consists of ultra hot dark matter. In the framework of structure formation, hot dark matter 
is disfavoured since in this case the galaxy-size density huctuations would get washed out 
by free-streaming leading to an earlier superclusters formation, while observations indicate 
early galaxy formation. 


Experiments involving dark matter particles detections, e.g DAMA 3^, search for parti- 
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cles with mass in the range (~ 15 - 120 GeV). More recently 3^, an upper bound limit for 
dark matter mass of ~ 197 TeV was established, based on relic abundance of thermal dark 
matter particles annihilating via a long-range interaction. 

Studies of the density perturbation evolution should be made in order to obtain an upper 
limit on the scalar field dark matter mass. We may perform this analysis in a future work. 
Another forthcoming study is trying to test this result against different potential dependen¬ 
cies of the scalar field, including self-interaction which could eventually furnish an upper 
limit to dark matter mass. 


VI. APPENDIX - THE MOTION EQUATION 


Another way to investigate the SFDM model is through the held motion equation. While 
it may be numerically intensive, this approach can be useful for dynamical interpretations. 
The motion equation for a scalar held is well known and can be obtained from the continuity 
equation as flT^ : 

rIV 

(45) 


• dV 

+ 3770 + — = 0 . 
dcp 


Writing it in terms of the e-folds number, N = Ina, which is more suitable for numerical 
integration, we have: 


+ 


77' 


^' + ^ = 0 , 


(46) 


where the primes at 0 and 77 denotes derivation with respect to N. Using the fact that 
can be written as P(f,{N) = i77^0'(iV)^ -f U(0), we can write the Hubble parameter as 


77" = 


SttG [U(0) pwe + Pa] 


(47) 


3 - 47rG0'(iV)2 

where pa =constant. Inserting this result into the motion equation, we arrive at an equation 
involving only held derivatives. However, in order to simplify the equations, we apply the 
free scalar held potential (U(0) = |m"0") and dehne the dimensionless quantities 


$ = 


SttG 


m 

^ ^ Ifo' 

With these dehnitions, the Friedmann equation can now be written 

^Ho) ~ 1 - i<F'(iV)2 


‘AO 


With this expression we may hnd the motion equation as 

1 


$" + 


2U2 


+ eJiAo) 4' + 2fi24] = 0, 


(48) 

(49) 

(50) 

(61) 
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where E‘^ = 


is given by fISOl) . This equation can be solved to hnd the free scalar held 


evolution. Furthermore, we have the dark matter density parameter 




3H^ 


_ $'2 


$2 


Thus: 


^(pO 


4 ? 

2 2 


(62) 


(53) 
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